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A  SYNERGISTIC  DAMAGE  MECHANICS  APPROACH  TO 
DURABILITY  OF  COMPOSITE  STRUCTURES 

1.  Goals  and  Objectives  of  Research 

The  overall  technical  goal  of  the  research  under  this  contract  has  been  to 
develop  a  methodology  for  assessment  of  durability  of  composite  structures 
operating  under  extreme  environments  typical  of  Army  vehicles.  The 
scientific  goal  has  been  to  combine  the  conventional  micromechanics  and 
continuum  damage  mechanics  approaches  into  a  synergistic  appoach  that 
maximizes  advantages  of  the  two  approaches  while  minimizing  their 
inherent  deficiencies. 

The  specific  objectives  set  for  the  three-year  program  have  been  to  apply  the 
concept  of  synergistic  damage  mechanics  to  composites  of  different 
configurations  in  order  to  validate  the  concept  as  well  as  to  help  formulate 
the  structure  of  the  methodology  for  life  assessment. 

The  proposed  research  was  planned  to  have  three  specific  tests.  Task  1 
dealing  with  mechanisms  of  damage  that  will  provide  the  physical  basis  of 
the  modeling  effort;  Task  2  focusing  on  formulation  of  stiffness-damage 
relationships;  and  Task  3  treating  the  evolution  of  damage.  These  tasks  have 
essentially  been  completed  and  the  objectives  set  have  been  satisfactorily 

met. 

2.  Summary  of  the  Most  Important  Results 
2.1  Accomplishments  of  Task  1. 

The  first  year's  research  in  the  program  focused  on  conducting  Task  1 
dealing  with  damage  mechanisms. 

The  first  laminate  configuration  taken  as  an  object  of  study  was  glass/epoxy 
[±6/902]s-  This  laminate  configuration  was  also  the  object  of  a  pilot  study 
that  was  used  as  a  motivation  for  the  proposal.  The  preliminary  results 
obtained  with  this  class  of  laminate  were  described  in  the  proposal.  Further 
results  were  obtained  during  the  course  of  this  project.  Those  results  were 
put  together  in  a  paper,  "An  Experimental  Study  of  Crack  Opening 
Displacement  in  Composite  Laminates  with  Varying  Constraints",  published 
in  International  Journal  of  Damage  Mechanics. 
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The  laminate  configuration  described  above  was  instrumental  in  validating 
the  concept  of  utilizing  the  crack  opening  displacement  (COD)  as  a 
micromechanics  parameter  in  the  second  order  tensor  for  characterization  of 
damage,  which  then  enters  as  the  internal  state  variable  in  the  continuum 
damage  formulation  of  material  response.  The  damage  occurring  in  the 
[±0/9O2]s  laminates  consists  of  transverse  cracks  in  the  90-plies  while  the  ±0- 
plies  provide  constraint  to  the  opening  of  those  cracks.  However,  the  off-axis 
plies  can  also  crack,  providing  another  mode  of  damage.  That  damage  can 
be  isolated  and  systematically  investigated  by  choosing  a  laminate 
configuration  such  as  [0,  ±04,  Oids-  This  was  done  and  the  results  were 
reported  in  the  paper,  "Damage  in  Composite  Laminates  with  Off-axis 
Plies",  published  in  Composites  Science  and  Technology. 

A  summary  of  the  results  reported  in  the  two  papers  mentioned  above 
follows. 

1.  "An  Experimental  Study  of  Crack  Opening  Displacement  in  Composite 
Laminates  with  Varying  Constraints" 

This  paper  reports  an  experimental  study  of  the  opening  displacement 
behavior  of  transverse  cracks  in  a  class  of  composite  laminates  where  the 
constraint  to  cracking  is  varied  systematically.  The  study  is  motivated  by  a 
synergistic  approach  to  damage  mechanics  [1]  where  the  crack  opening 
displacement  forms  a  key  micromechanically  determined  parameter  in  a 
continuum  damage  formulation  of  the  composite  response.  It  is  found  that 
although  the  opening  displacement  shows  scatter  from  point  to  point  along  a 
crack  length,  its  average  over  the  crack  length  forms  a  robust  parameter. 
This  parameter,  when  introduced  as  the  damage  influence  coefficient  in  a 
continuum  description  of  damage,  allows  satisfactory  predictions  of  the 
elastic  moduli  changes  due  to  transverse  cracks  in  the  considered  class  of 
laminates  with  varying  constraints. 

Highlights  of  Results: 

Figure  1  shows  the  normalized  average  value  of  COD  of  transverse  cracks 
plotted  against  0  in  the  [±0/9O2]s  laminates.  It  is  apparent  that  the  constraint 
to  COD  is  highest  at  0  =  0,  i.e.,  for  cross  ply  laminates,  and  it  diminishes 
with  increasing  0,  vanishing  at  0  =  90.  This  COD  provides  a  link  between 
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stiffness  changes  due  to  transverse  cracking  in  [±0/9O2]s  laminates  at 
different  values  of  0.  Defining  a  constraint  parameter  in  the  CDM 
formulation  of  stiffness-damage  relationships,  and  setting  it  in  proportion  to 
the  normalized  COD,  allows  prediction  of  stiffness  changes  due  to 
transverse  cracking  in  [i0/9O2]s  laminates  when  the  stiffness  changes  in  the 
cross  ply  laminate  are  known.  The  predictions  by  this  procedure  are 
compared  with  experimental  data  in  Figs.  2-4. 


Fig.  1.  Variation  of  the  normalized  COD  with  6  in  [±0/9O2]s  laminates. 
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Fig.  2.  Predictions  using  experimental  COD  values  of  (a)  longitudinal  Young's  modulus 
and  (b)  Poisson's  ratio  compared  with  experimental  data  for  [+15/902]$  laminates. 
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Fig.  3.  Predictions  using  experimental  COD  values  of  (a)  longitudinal  Young's  modulus 
and  (b)  Poisson's  ratio  compared  with  experimental  data  for  [±30/902]s  laminates. 
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Fig.  4.  Predictions  using  experimental  COD  values  of  (a)  longitudinal  Young's  modulus 
and  (b)  Poisson's  ratio  compared  with  experimental  data  for  [±40/902]s  laminates. 
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Conclusions:  The  applicability  of  continuum  damage  mechanics  to  a  wide 
class  of  laminate  configurations  and  crack  geometries,  as  against 
micromechanics,  which  is  laminate  and  crack  geometry  specific,  is  at  the 
cost  of  the  necessity  to  determine  material  constants  experimentally.  This 
situation  can  be  improved  by  a  synergistic  damage  mechanics  approach 
where  judiciously  selected  micromechanics,  performed  in  a  limited  measure, 
provides  input  into  a  continuum  damage  mechanics  formulation. 

For  the  class  of  laminates  given  by  [±0/9O2]s,  we  have  demonstrated  here 
that  the  average  crack  opening  displacement  of  90-ply  cracks  provides  the 
micromechanics  parameter  which  relates  well  with  the  material  constant 
which  varies  with  the  constraining  ply  angle  0.  Thus,  use  of  this 
micromechanics  parameter  will  avoid  the  necessity  to  experimentally 
determine  that  material  constant,  which  usually  is  a  demanding  task 
involving  testing  of  a  large  number  of  specimens  for  each  value  of  0. 

The  crack  opening  displacements  of  90-ply  cracks  were  determined  here  by 
actual  physical  measurements.  This  was  done  to  gain  insight  into  the 
phenomenon  of  constraint  effects  on  crack  opening  displacement,  in  addition 
to  validating  it  as  a  micromechanics  parameter  for  the  purpose  stated  above. 
It  is  obviously  not  always  feasible  to  determine  this  parameter  by  physical 
measurements.  As  we  have  suggested  previously,  this  parameter  should 
instead  be  computed  by  an  approximate  numerical  scheme  such  as  the  finite 
element  method.  This  was  accomplished  in  a  later  work  described  below. 

Finally,  it  is  noted  that  the  treatment  presented  here  has  been  for  the  case  of 
a  linear  dependence  of  moduli  changes  with  crack  density.  For  the  cases 
where  nonlinearity  is  displayed,  as  crack  densities  approach  their  saturation 
values,  higher  order  terms  must  be  included  in  the  stiffness-damage 
relationships  [2]. 


2.  "Damage  in  Composite  Laminates  with  Off-Axis  Plies" 

Damage  in  off-axis  plies  of  composite  laminates  is  studied  by  examining  the 
configuration  [0,  +04,  0i/2]s  with  0  =  25,  40,  55,  70  and  90  subjected  to 
tensile  loading  in  the  axial  direction.  It  is  found  that  for  the  values  of  0 
where  the  stress  in  the  off-axis  plies  normal  to  the  fibers  is  tensile,  ply  cracks 
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lying  along  fibers  initiate  and  increase  in  number,  while  for  other  0  values 
the  plies  do  not  undergo  this  damage,  as  expected.  However,  the  overall 
laminate  elastic  moduli  are  found  to  change  also  for  the  0  values  where  no 
ply  cracks  exist.  It  is  postulated  that  a  shear-induced  degradation  of  the  off- 
axis  plies  is  responsible  for  the  observed  laminate  moduli  changes.  The 
prediction  of  changes  in  these  moduli  by  using  the  ply  shear  modulus 
measured  on  [+04]s  appears  to  support  this  postulate.  For  the  case  of  moduli 
changes  caused  by  ply  cracks  the  recently  proposed  synergistic  damage 
mechanics  approach  [1]  is  applied.  The  implications  of  the  findings  of  this 
work  on  a  class  of  continuum  damage  mechanics  formulations  proposed  in 
the  literature  are  discussed. 

Highlights  of  Results: 

The  average  CODs  of  [0,  ±04,  Oids  laminates  for  0  decreasing  from  90  are 
shown  in  Fig.  5.  The  interesting  result  is  that  the  COD  vanishes  at  0  <  40. 
For  these  cases  the  stress  normal  to  fibers  in  the  0-plies  is  compressive. 
However,  the  laminate  stiffness  properties  are  found  to  change  whether  the 
off-axis  plies  crack  or  not.  When  cracking  occurs,  the  moduli  changes  are 
predicted  by  the  synergistic  damage  mechanics  procedure  described  in  the 
previous  paper.  Predictions  thus  made  are  compared  with  experimental  data 
in  Fig.  6-7  for  [0,  ±704,  0i/2]s  laminate.  For  [0,  ±554,  0i/2]s  laminate  it  is  found 
that  the  ply  cracks  are  only  partly  responsible  for  the  laminate  stiffness 
changes.  Figures  8-9  shows  how  the  stiffness  predictions  due  to  cracking, 
marked  as  CDM,  do  not  agree  with  the  experimental  values.  The  other 
effect,  responsible  for  the  remaining  stiffness  changes,  is  due  to  the  shear 
stress  induced  damage.  Predictions  due  to  that  effect  as  well  as  the  total 
effect  by  addition  of  the  two  are  also  shown  in  the  figure.  Finally,  Figs.  10- 
1 1  compare  the  stiffness  changes  due  to  the  shear  stress  induced 
degradation,  when  no  ply  cracks  are  present,  for  [0,  ±404,  0i/2]s  laminate. 

Conclusions 

Based  on  a  systematic  investigation  of  damage  and  the  associated 
deformational  response  of  laminates  of  configuration  [0,  ±04,  0i/2]s  with  0  = 
25,  40,  55,  70  and  90  loaded  in  axial  tension  it  has  been  found  that  two 
different  damage  modes  exist  depending  on  the  angle  0.  For  the  case  where 
the  ply  stress  component  normal  to  the  fibers  is  tensile,  the  familiar 
mechanism  of  multiple  matrix  cracks  along  fibers  exists  and  the  laminate 
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elastic  moduli  change  as  a  function  of  the  density  of  these  cracks.  A 
procedure  for  predicting  moduli  changes  based  on  the  synergistic  damage 
mechanics  proposed  earlier  [1]  is  found  to  be  successful  for  this  case.  When 
the  stress  component  normal  to  fibers  is  compressive  the  multiple  matrix 
cracking  does  not  occur,  but  the  elastic  moduli  still  undergo  changes.  It  has 
been  demonstrated  that  the  underlying  cause  of  these  moduli  changes  is 
damage  induced  by  shear  stresses  in  off-axis  plies.  A  procedure  for 
predicting  these  moduli  changes  has  been  proposed  and  demonstrated  to 
work  well. 

The  nature  of  the  shear  stress  induced  ply  damage  is  a  subject  of  further 
research.  Its  existence  raises  questions  concerning  applicability  of  certain 
continuum  damage  models  that  define  damage  in  terms  of  the  moduli 
changes. 


Fig.  5.  The  average  COD  plotted  against  the  orientation  angle  9  in  [O/±04/O>/,]s  laminates 
at  an  applied  axial  strain  of  0.5%. 
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Fig.  6.  The  longitudinal  Young’s  modulus  normalized  by  its  virgin  state  value  plotted 
against  the  crack  density  in  [0/±704/0i/2]s  laminate.  The  solid  line  is  the  prediction  of  the 
CDM  model. 


Fig.  7.  The  Poisson’s  ratio  normalized  by  its  virgin  state  value  plotted  against  the  crack 
density  in  [0/±704/0i/2]s  laminate.  The  solid  line  is  the  prediction  of  the  CDM  model 
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Fig.  8.  The  longitudinal  Young’s  modulus  normalized  by  its  virgin  state  value  plotted 
against  the  crack  density  in  [0/±554/0i/2]s  laminate.  The  dotted  line  is  the  prediction  of  the 
CDM  model. 
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Fig.  9.  The  Poisson’s  ratio  normalized  by  its  virgin  state  value  plotted  against  the  crack 
density  in  [0/±554/0i/2]s  laminate.  The  dotted  line  is  the  prediction  of  the  CDM  model 
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Fig.  10.  The  longitudinal  Young’s  modulus  normalized  by  its  virgin  state  value  plotted 
against  the  applied  axial  strain  in  [0/±404/0i/2]s  laminate. 
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Fig.  11.  The  Poisson’s  ratio  normalized  by  its  virgin  state  value  plotted  against  the 
applied  axial  strain  in  [0l±A0\IQ\n\s  laminate. 
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2Jl  Accomplishments  of  Task  2 

The  research  during  the  first  year  (6/1/97  -  5/31/98)  was  focused  on  Task  1 
dealing  with  damage  mechanisms  caused  by  mechanical  loading  ot 
composite  materials.  Two  configurations  of  composite  laminates  were 
considesred,  namely,  [±0/902],  and  [O/±04/O%],  with  0  *=  25, 40,  55,  70 
degrees.  These  configurations  were  selected  to  reveal  and  to  chai  acte^  me 
constraint  effects  on  crack  surfeoe  displacements  and  the  nature  of  damage 
in  off-axis  plies.  These  objectives  were  accomplished  during  the  first  year  s 
work.  In  order  to  continue  on  to  Task  2  aimed  at  formulating  stifmess- 
damage  relationships  in  a  synergistic  damage  mechanics  approach  we 
needed  to  develop  micromechanics  analyses  of  crack  surface  displacements 
as  well  as  evaluate  the  shear  stress  induced  off-axis  ply  damage  ^at  was 
discovered  during  the  first  year’s  research.  This  was  accomplished  as 

described  below. 


A  systematic  computational  study  was  conducted  to  study  the  effects  of  the 
constraming  (±6)  plies  on  crack  opening  displacement  (COD)  of 
cracks  in  laminates.  Replacing  the  ±0  plies  by  an  equivalent  sublaminate  (S) 
the  COD  behavior  in  [S^9Q^,  laminates  was  staged  by  varying  p^etws 
such  as  the  sublaminate  stififeess  and  thickness  (i.e.  numbers  n  and  m)  o 
the  90-plics.  The  three  oases  of  plane  stress,  plane  strain  and  general^ed 
plane  strain  were  considered.  It  turned  out  that  they  all  produced  essentially 
the  same  COD.  Figures  12-14  show  the  variation  of  COD  with  the 
stiffiiess  (longitudinal  Young's  modulus)  of  the  sublaminate  S.  The  COD 
plotted  in  these  figures  is  normalized  by  the  lon^tudinal  strain  and  the  half 

thickness  (d)  of  90-plie3. 


Pig.  12.  The  calculated  variation  of  COD  of  90-pJy  ctacks  in  [S,.  904]*  laminates  >with  the  wWivn  «iffhess 
of  sublaminaxc  S. 
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Fig.  13.  The  calculated  variation  of  COD  of  90-ply  cracks  in  [S4,  904]s  laminates  with  the 
relative  stiffness  of  sublaminate  S. 


Fig.  14.  The  calculated  variation  of  COD  of  90-ply  cracks  in  [Sg,  904]s  laminates  with  the 
relative  stiffness  of  sublaminate  S. 
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Comparison  of  these  three  figures  provide  several  interesting  results: 

1.  For  highly  anisotropic  composites  such  as  carbon-epoxy  cross  ply 
laminates,  where  the  ratio  of  the  0-ply  to  90-ply  stiffness  is  large,  the 
constraining  ply  thickness  plays  a  minor  role  in  affecting  the  COD  of  the 
90-ply  cracks  normalized  by  the  90-ply  thickness. 

2.  For  constraining  ply  orientation  0  deviating  significantly  from  the  axial 
direction,  say  by  20  degrees  or  more,  and  for  moderately  anisotropic 
laminates,  e.g.,  glass-epoxy  cross  ply  laminates,  the  constraining  ply 
thickness  has  a  significant  effect  on  COD  of  90-ply  cracks. 

3.  For  all  cases  studied  the  normalized  COD  tends  asymptotically  to  1.0  for 
large  normalized  sublaminate  to  90-ply  stiffness  ratio. 

It  is  to  be  noted  that  the  results  presented  here  are  for  non-interacting  cracks. 
Previous  studies  by  this  author  and  his  associates  [3,4]  have  indicated  that  in 
most  cases  linearized  stiffness-damage  relationships  suffice,  in  which 
instance  CODs  of  non-interacting  cracks  will  provide  the  results  needed. 

The  expressions  shown  in  Figs.  12-14  provide  functional  form  for  variation 
of  COD  that  can  be  incorporated  in  this  author's  synergistic  damage 
mechanics  formulation  [1].  Preliminary  comparisons  of  stiffness  changes  by 
use  of  these  expressions  and  by  use  of  experimentally  determined  CODs  [4] 
show  better  agreement  with  experimental  data.  A  paper  is  under  preparation 
to  report  this  comparison. 

2.3  Accomplishments  of  Task  3 

The  work  in  Task  2  had  focused  on  determining  the  changes  in  the  stiffiiess 
properties  at  fixed  states  of  damage.  The  work  in  this  task  combines  this 
with  prediction  of  damage  evolution.  Thus  the  stiffness  properties  are 
described  as  functions  of  the  applied  loading.  The  specific  case  considered  is 
that  of  transverse  ply  cracking  constrained  by  sublaminates  of  varying 
stiffness.  Figure  15  shows  the  laminate  with  cracks  in  the  90  degree  layers 
sandwiched  between  sublaminates  of  ±0  plies.  The  angle  0  is  varied  to  have 
values  0,  15,  30  and  40  degrees.  The  previously  described  constraint 
parameter  is  calculated  by  a  parametric  study  in  which  the  ply  thickness  and 
ply  properties  are  changed.  A  relationship  is  developed  to  describe  the 
effects  of  relative  ply  stiffness  and  the  ratio  of  crack  size  to  the  constraining 
ply  thickness.  A  fracture  mechanics  based  procedure  is  developed  to  predict 
the  evolution  of  ply  cracking.  This  procedure  involves  Monte-Carlo 
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simulation  incorporating  the  statistical  variation  of  the  fracture  toughness  in 
the  fiber  direction. 

Figures  16  and  17  illustrate  the  predictions  of  deformational  response  (in 
tenns  of  the  axial  Young’s  modulus  and  Poisson’s  ratio)  with  applied  strain 
for  the  cases  of  0  =  15  and  0  =  40. 


Fig.  15.  A  composite  laminate  loaded  in  tension  axially,  developing  multiple  cracks  in  the 
middle  90  degree  layers  with  the  outer  sublaminate  S  remaining  uncracked. 
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Fig.  17.  The  variations  of  the  axial  Young’s  modulus  and  Poisson’s  ratio  with  the  applied 

strain  for  0  =  40. 
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